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We present results on an analysis of the decay constants fsd and /g, with two flavours of sea quark. The 
calculation has been carried out on 3 different bare gauge couplings and 4 sea quark masses at each gauge 
coupling, with m^/m p ranging from 0.8 to 0.6. We employ the Fermilab formalism to perform calculations with 
heavy quarks whose mass is in the range of the b-quark. A detailed comparison with a quenched calculation using 
the same action is made to elucidate the effects due to the sea quarks. 



1. Introduction 

An accurate determination of the parameters 
fBd\/B B d and £ = f B W B Ba / fBdV B Bd , in con- 
junction with the (future) experimental data 
on Amd (and Am s ) will provide excellent con- 
straints on \Vtd\ an d |Vt a |/|Vtd|(l). Their calcu- 
lation in the quenched approximation, using the 
plaquette gluon action, has been carried out us- 
ing a number of different formulations for heavy 
quarks and the results are converging 

A major uncertainty in these results, however, 
is that they may be susceptible to large correc- 
tions due to the effect of sea quarks @|e|. The 
penultimate step in eliminating this uncertainty is 
to consider the effect of two flavours of sea quark 
(Nf = 2). Here we present the status of such 
a study, using the clover action for heavy quark, 
which we started last year||. 

As the lattice spacings available to us are rel- 
atively coarse, we deal with the large mass of b 
quark rriba > 1 in the formalism of Ref. [Q, which 
has been previously applied in the quenched cal- 
culation of B meson decay constants pjj^j. 

In order to reduce discretisation effects, we 



Table 1 

Sizes of lattices used in calculation. 



Y 3 x N t 


P 


■*V s&chiral 


^chiral 


12 3 x 24 


1.8 


2.8 fm 


0.835 GeV 


16 3 x 32 


1.95 


2.5 fm 


1.263 GeV 


24 3 x 48 


2.1 


2.6 fm 


1.807 GeV 



*talk presented by H.P. Shanahan 



have employed an RG-improved action in the 
gluon sector. Since this action was not considered 
in previous studies of fs , we repeat the calcula- 
tion in quenched QCD to compare with full QCD 
results. 

A comparison is also made with preliminary 
NRQCD results for the decay constants obtained 
on the same full QCD configurations, presented 
elsewhere in these proceedings [fTof . 

2. Computational Details 

Gauge configurations are computed using 
the RG-improved gauge action and a tadpole- 
improved SW clover fermion action to create the 
background for two flavours of sea-quark. As 
listed in Table |l|, three bare gauge couplings were 
used. At each bare gauge coupling (/3) runs were 
made for four sea quark masses {m s q ea ) in the 
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range m^/rrip w 0.8 — 0.6. More details of these 
gauge configurations are presented in Ref . • 

At each ((3,m s q ea ), propagators for 8 heavy 
and 2 light valence quark masses were computed. 
The same tadpole-improved clover action is em- 
ployed for both heavy and light quarks. Details 
of the propagators computation and the choice of 
smearing can be found in Ref. [||. 

The heavy quark masses ranged roughly from 
charm to bottom. One of the light valence quark 
masses was chosen so that mps/rnv = 0.688, 
which approximately corresponds to the strange 
quark mass. The other light valence quark mass 
was chosen equal to the bare sea quark mass pa- 
rameter. For each (3, both types of matrix ele- 
ment were fitted as a function of (m* ea ) 2 to com- 
pute fg s and fsd at zero sea quark mass and finite 
lattice spacing. 

We include the correction to the axial vector 
current obtained by the following substitution to 
the heavy quark field , 



h{x) ->■(!- diaj.D)h(x) 



(1) 



Using tadpole-improved tree- level value for d\ , we 
find that this correction has a contribution of up 
to 8% of the original current. 

As a definition for the ground state mass of 
the heavy-light meson we used the "HQET" mass 
definition, as discussed in Ref. []l3|,PU| . The phys- 
ical scale of lattice spacing is fixed by the string 
tension measured for each sea quark mass using 
V^= 427 MeV. 

We made the quenched calculation with the 
same action at 5 different gauge couplings each 
on 24 3 x 48 and 16 3 x 32 lattices. The couplings 
were tuned so that they match the string tensions 
of full QCD for the four values of finite sea quark 
mass and the extrapolated zero sea quark mass. 
For more details see Ref. . 

3. Results 

In Fig. [l] we plot results for fs s for Nf = 2 
and in the quenched approximation (Nf — 0); 
sea quark mass is extrapolated to the chiral limit 
for Nf = 2. A clear increase of 10-20% is seen 
from two flavours of sea quark. 

We also include the preliminary results from 
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Figure 1. A comparison of fg s for full and 
quenched QCD. NRQCD results are also shown 
(squares). The scale is set by string tension 
yfd = 427 MeV. 
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Figure 2. A comparison of preliminary NRQCD 
and clover results for 

NRQCD |§ in Fig. @, for two values of finite 
sea quark mass at (3 = 1.95 in full QCD (filled 
squares), and one value in quenched QCD (open 
squares). In Fig. ^ we present a similar compari- 
son of the Fermilab and NRQCD approaches for 
fsd- In both figures we find good consistency of 
results between the two approaches. 

Finally, we plot the ratio /s 3 //s d in Fig. || us- 
ing the K meson mass to set the strange quark 
mass. We observe only mild variation of the ra- 
tio with respect to a for even the coarsest of our 
lattice spacings. 

As a preliminary result we quote 



J E 



= 251 ± 3 ± 4(m s ) ±$ (fit) MeV, (2) 
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• n,=2 Clover 

■ n,=2 NRQCD, p=1 .95 
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Figure 3. The fs„/ fs d as a function of a 2 a for 
the relativistic and NRQCD actions. 



f n f- 

J B 

f n f-~ 
J Bs 



f n f =2 
JB 



210 ± 7 ±\ A (fit) MeV, (3) 
1.20 ± 4 ± 2(m s ) ±| (fit) . (4) 



The first errors are statistical. The error labelled 
(m s ) is due to the ambiguity of using the mass of 
the (j) or K to set the strange quark mass. The 
central values are determined by assuming the re- 
sults independent of a for the two finer lattice 
spacings. The resulting systematic error (fit) is 
derived by taking the difference between a con- 
stant and a linear fit in a for all three points. We 
should also add an uncertainty due to the choice 
of scale; our preliminary estimate is 15-20 MeV 
from comparision of results using the scale deter- 
mind from m p . 

4. Conclusions 

At this point it seems clear that there exists 
a systematic difference between the Nf — and 
2 data. Encouragingly enough, the preliminary 
NRQCD results are also in agreement with the 
relativistic results in both cases as well. One wor- 
rying point is that the quenched results for fs s 
is approximately 10% larger than the quenched 
results using the Wilson action (|||. Clearly 
this effect needs further examination. The ratio 
fBs/ fBd appears to be less affected by discreti- 
sation effects and is not substantially different 
from previous quenched calculations. It should be 
noted that even a systematic error of 10% for this 
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ratio would be of substantial use for phenomenol- 
ogists. It seems plausible then that a calculation 
of fBs/fBd could be carried out for three flavors 
of dynamical quarks on a comparatively coarse 
lattice. 
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